Rayleigh- Taylor turbulence is nothing like Kolmogorov's in the self similar regime 
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An increasing number of numerical simulations and experiments describing the turbulent spectrum 
of Rayleigh- Taylor (RT) mixing layers came to light over the past few years. Results reported in 
recent studies allow to rule out a turbulence a la Kolmogorov as a mechanism acting on a self similar 
RT turbulent mixing layer. A different mechanism is presented, which complies with both numerical 
and experimental results and relates RT flow to other buoyant flows. 
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A Raylcigh-Taylor (RT) instability occurs when- 
ever a light fluid; pi, pushes a heavy fluid, p2, or simi- 
larly, when a heavy fluid on top of a lighter fluid is sub- 
ject to a gravitational acceleration. The understanding 
of such instability in the developed turbulent regime is 
of primary interest to many fields of physics and tech- 
nology since it is an important cause of mixing between 
two fluids of different density. In astrophysics for in- 
stance, it is responsible for the outward acceleration of 
a thermonuclear flame in type la supernovae Q, but it 
also plays an important role in shaping the interstellar 
medium so that new stars can be born [J| . The technol- 
ogy of confinement fusion also relies on a good under- 
standing of RT mixing [Bl and ways to reduce it Q . The 
RT flow of two incompressible fluids in the low Atwood 
limit, A = {p2 — Pi)/{p2 + Pi) ^ 1 (Boussinesq approx- 
imation), is governed by a concentration equation (fT|). 
the Navier Stokes equation supplemented with a buoy- 
ant source term ^ and the incompressibility constraint 
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where g is a stationary and uniform gravitational acceler- 
ation vector fleld (i.e. planar symmetry is assumed). The 
coefficient k is the molecular diffusion coefficient and ly 
is the kinematic viscosity of the mixture. They are both 
supposed constant. Without loss of generality, g is par- 
allel to the z-axis. This is why, for any generic physical 
value <&, the average (so defined numerically and experi- 
mentally) will be ($)(z) = Js ^ throughout this 
paper. The fluctuating part will be denoted with a prime 
and defined as = $ — ($). 

With the increasing capacity of super computers, 
many simulations of RT flows in the developed turbu- 
lent regime have been performed, which describe the ve- 
locity spectrum E{k) 0, S, Q, defined in such a way 
that (m'^) = / dkE{k), or the concentration spectrum 
i^c(A;)Jl0jn|, El, m, defined as (c'^) = / dkEc{k), or 
both [3]i|y In the same time, although fewer in num- 



ber, experimental investigations of E{k) [I3 and Ec{k) 
[lol . IitI . [lit have been carried out. A quick inspection 
of these results shows that no consensus arises concern- 
ing the shape of these spectrum. From a theoretical 
point of view, the situation is not more satisfactory. In 
[l9| it is claimed that the Kolmogorov-Obukhov scheme, 
E{k) ~ holds in 3D RT mixing given that the ef- 

fect of buoyancy on turbulence, although fundamental at 
the largest scale, becomes irrelevant at smaller scales. In 
[20t the particular RT time scale l/y'kgA at wave number 
k has been postulated to vary like the turbulent spectral 
transfer time scale e k~'^ E~'^{k) yielding an RT turbulent 
spectrum E{k) ~ k^"^!^. Another troublesome point is 
that there are no convincing objective criteria to assert 
whether or not an RT flow has reached the gravitational 
self similar regime but to plot the mixing zone width L(i) 
versus A g and sec if there is a straight line somewhere 
without any information on the expected slope. Subse- 
quently, suggesting the possibility of any behavior of the 
velocity spectrum in the self similar regime and compar- 
ing it to experiments and simulations can be dubious if 
it is not known for certain that these experiment or sim- 
ulations have reached this regime. So, the flrst aim of 
this paper will be to describe such an objective criterion. 
The ultimate goal will be to prove that the Kolmogorov 
mechanism does not explain the observed numerical and 
experimental results. A theory, based on a spectral equa- 
tion, will be presented which shows that a balance mech- 
anism between buoyancy and spectral energy transfer can 
settle at low wave numbers in the self similar regime. 

The RT developed turbulent regime, as complex as it 



may look, is thought to evolve self similarly [2l|. The 
size of the largest significant turbulent structure, there- 
fore, grows like the size of the turbulent mixing zone 
which evolves as L{t) — aAgt^ where a is the 
mixing-zone-width growth rate parameter whose value 
is around ~ 0.06 experimentally and between 0.03 and 
0.07 numerically. It can be assumed that at low wave 
numbers in the self similar regime the velocity spectrum 
E{k, t) ~Q for k < ki{t), where ki{t) is the wave number 
corresponding to the maximum of the velocity spectrum 
(this point will be called A due to the shape of the ideal- 
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ized spectrum at this location), and E{k,t) w ^i{t) fc~"' 
for k > ki(t). Nothing is assumed concerning the be- 
havior of E(k,t) at intermediate and high wave num- 
bers k ^ ki{t). In the self similar regime, fc;(<) de- 
creases like ~ {Agt^^ and the mean turbulent kinetic 
energy increases like (u'^) ^ {Agtf' [2l[. Provided that 
the spectrum decreases toward high wave numbers up to 
k = krj{t) above which it is zero, it can be found that 



dkE{k,t) 
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if ri; > 1 and kjj ^ ki. We could have taken into account 
an additional and more realistic ~ fc^ spectrum in the 
region k < ki (t) instead of zero but that would only have 
changed the unimportant constant coefficient in front of 
the above result. For the previously depicted self simi- 
lar evolution to occur, using (g]), the spectrum level at 
low wave number must verify ^ (-^ff)^ 
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The parametrized trajectory of A in the 



k-E plot can then be determined since kx{t) ~ ki{t) and 
Ex{t) — E{kx{t),t). After eliminating the variable t, it 
is found that the point A must evolve on the curve 
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independently of the slope of the spectrum (rij). This is 
a universal condition in the sense that both those who 
believe in a Kolmogorov scenario ni = 5/3 and others 
will agree on. It is objective since it amounts to look 
for a straight fine in the log(£^)-log(A:) plot whose slope, 
this time, is unambiguously determined (—2). Therefore, 
such behavior should be checked whenever an RT flow is 
said to have reach a self similar regime. This regime 
requires large resolution so that the mixing zone can ex- 
pand and reach the self similar regime before it collides 
with the simulation box border. Indeed, large resolution 
simulations [1, 14, 1^ show the mark ([5]) of the self sim- 



ilar regime which starts near the end. The only experi- 
mental result which shows the evolution of the velocity 
spectrum [l6l | also complies with this rule. 

It is now possible to select, among all the simulations 
cited in the introduction, those which have reached the 
self similar regime. If it is assumed that RT turbulence 
follows Kolmogorov's mechanism IIH, the turbulent 
spectrum has the well known form E(k) = Ck k~^/^ . 
Using Eq.Q with = 5/3 > 1 and the self similar laws 
for ki{t) and (u'^), it is straightforward to conclude 
that e ~ t and fc^(t) ~ [e/v'^^^^'^ ^ t^^^ and also that 
\I';(t) = Cif e(t)2/3 t^/^. This conclusion is refuted by 
recent DNS/LES simulations [l,[l3l showing the velocity 
spectrum evolution in time. Even though it is not stated 
in these references, it can clearly be noticed that the level 
of the velocity spectrum at low wave number, 4"; (t) , does 
not grow as time evolves but remains still. This is an 
important observation which constrained the mechanism 
of RT turbulence. 



In order to understand the mechanism, it is worth writ- 
ing the equation governing the evolution of E{k,t) (a 
generalisation of Lin's equation [i^]) out of the averaged 
second moment equation of ([2|) to make the buoyancy 
production term appear. It can be achieved by multi- 
plying (121) by M. The spectral equation can then be de- 
rived by applying a Fourier transform, ~ in the xy plan 
to the resulting equation and retaining the zero mode 
contribution. In addition to the rhs of Lin's equation, 
we find a spectral buoyancy production term deriving 
from 2^5<c'(0) = Ag J (Pk [u,{k)c* (k) +ul{k)c{k)] 
which, assuming homogeneity and isotropy in the xy 
plan and phase coherence between Uz{k) and c(fc), 
reduces to / dk Ag eI^^ {k)E^/'^{k) since ~ 
k-y^E^/^k)/V2^ and \c{k)\ - k"^/^ eI^'^ {k)/V2^. 
The resulting generalisation of Lin's spectral equation is 
therefore 



dtE = T{k) - 2vk^E -f pAgEl'^E^''^ 



(6) 



where (3 accounts for phase incoherence between c and 
Uz ■ It depends on k but remains of order unity. The first 
term in the rhs (T term) is the so called spectral energy 
transfer. In the case of a forced turbulence (Kolmogorov 
mechanism), this term is negative in the low wave num- 
ber, to balance with the production, it is approximately 
zero throughout the inertial range and become positive in 
the dissipative range to balance with dissipation. It ac- 
counts for the non linear triad interaction responsible for 
the forward cascade in forced turbulence and must verify 
dkT{k) — 0. The second term is the exact contribu- 
tion of viscosity in the spectral evolution equation and 
it is responsible for dissipation at high wave number in 
forced turbulence. The last term, on the other hand, 
has never been written to the author's knowledge. It is 
this term which makes RT flows and buoyant flows in 
general so different. As expected intuitively, it depends 
on the concentration spectrum because it is a concen- 
tration heterogeneity that induces motion through buoy- 
ancy. In order to distinguish the influence of all three 
terms, knowledge of the order of magnitude of T(fc) is 
required. This term has contribution from u'ju'^ which 
in spectral language means E^/"^ (power counting of u' 
is 3 which means 3/2 in term of E) and more precisely 
/j3/2^3/2 homogeneity reason. The same sort of ar- 
gument is used with the pressure term u'^p' which is a 
non local term with two sources : pure advection and 
buoyancy. Thus, the pressure contribution brings an- 
other k'^l'^E^^'^ due to pure advection and a AgE^^ E^^'^ 
for buoyancy (power counting of u' and c' in buoyancy 
yields 1/2 for E and Ec). And flnally, vdju'^ is responsi- 
ble for a vk^E contribution. In a nutshell, T{k) has three 
spectral contributions whose order of magnitude are (i) 
A:3/2£;3/2^ (n) vk'^E and (in) A gE^^E^/^. Thus, the evo- 
lution of the velocity spectrum density in Eq.Q is also 
controlled by these three contributions. It is enlighten- 
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FIG. 1: Predominance diagram of various terms in tlie spec- 
tral equation ((6]). In tliis diagram, tiie point A evolves on the 
A line following ki{t). The gray curve represents the velocity 
spectrum in the RT self similar regime. 



ing to draw on a plot (Fig[T]) the predominance domains 
of the non linearity (i), viscosity (ii) and buoyancy (iii) 
with k on the horizontal axis and E(k) on the vertical 
axis. It is done by equating the three terms two by two. 
It gives three boundary lines in a log-log plot which will 
be referred to as the non-linearity-buoyancy (NLg), the 
non-lincarity-viscosity (NL^) and the buoyancy-viscosity 
(g:^) boarders respectively described by 

^NLg(fc,t) ~ Agk-'^/^E.ik^t)'/^ , (7) 
E^i^^ik) ^ j^H , (8) 
E^^{k,t) - {Agf 1^^^ k-* Ecik,t) . (9) 

Firstly, Ec{k,t) may not have the same power law at 
low wave number (NLg) in ([7]) and at high wave num- 
ber (NLi/) in ([9]). In the self similar regime (c' ) tends 
to be a constant [2l| which will be denoted Cn. At 



c'^), it is 



low wave numbers, by applying Eq.(|4]) to 
found that Ec{k,t) must vary like clki"-\t)k-"^^ ~ 
cKAg)^^""^ t'^i^-"'") k'"^"' , where rid is undetermined 
yet. That is why, in the self similar regime and at low 
wave number, Eq. ([7]) may be refined as 



Enhgik) - Co {Ag) ^ t 



(10) 



It is to be noticed that when rid = 1; tl^c previous time 
evolution is changed to \og{t). Secondly, it is worth notic- 
ing that £'NLi/(fc) ^ k is independant of the concentra- 
tion spectrum Ec{k,t), that it docs not depend on time 
and that it is the exact line where the inertial cascade 
vanishes and where dissipation starts acting in the Kol- 
mogorov mechanism. 

It is now possible to describe the evolution of the ve- 
locity spectrum in an RT mixing layer. Let us imagine 



that a peaked initial condition is chosen at a wave num- 

1 /3 

ber below (^Ag/v^) . Spontaneous RT can then oc- 

1/3 

cur. The closest to (^Ag/v'^) , the fastest is the linear 
growth. The spectrum grows until it reaches the NL do- 
main where mode coupling can start. At the beginning, 
since both E{k) and Ec{k) are peaked around the initial 
wave number, so is the buoyancy production which then 
acts like a narrow band forcing. A Kolmogorov spec- 
trum can then settle between the initial wave number 
and Nhv (where energy is dissipated). This is precisely 
what is observed after the linear growth regime in the 
mixing phase. As time goes by, both spectrums spread 
over a wider range of wave numbers and when mixing is 
established, production becomes broad band and spec- 
tral energy transfer must balance with buoyancy at the 
lowest wave number. That is why, the velocity spectrum 
lies alongside i?NLg(^) (sec FiglT]), given by ([7]) and also 
(fTO|l in the self similar RT case. 

This balance mechanism is observed for Rayleigh- 
Bcnard (RB) flows for it explains its velocity and con- 
centration spectrum layout. RT and RB flows are both 
governed by the same set of equation (concentration is 
replaced by temperature in RB) but boundary condi- 
tions in RB are independent of time and so are the veloc- 
ity and concentration spectrum. At low wave numbers, 

dtc''^ ^ ~dj ^u^c'^^ which in spectral language approx- 
imately means dtEc ~ dk {k^/"^ E^/"^ Ec) ■ If Ec is inde- 
pendent of time at low wave numbers in the self similar 
regime, that means the spectral transfer of concentration 
'--^ k^/^E^/^Ec must be independent of k, that is to say 
| — ric/ = for RB flows (power counting of /c). More- 
over, since (O must be valid in RB flows (Eq. ((TO)) is not 
valid in this case because RT boundary conditions have 
been used to derive it) the relation — n; = — | — ^ must 
also be true. Together with the previous relation it is 
then foimd that n; = = 2.2 and rid = | = 1-4, which 
is exactly the Bolg iano-Obukov (BO) [1^ scaling found 
in RB flows 2J,|25|. This mechanism is also corroborated 
by experiment [l6| where velocity spectrum level at low 
wave numbers can be seen to decrease in time although 
the turbulent kinetic energy increases. That means, pO[) . 
that rid ^ li which is confirmed by all numerical and 
experimental results, and that, once again, Kolmogorov 
mechanism is ruled out since it predicts a velocity spec- 
trum level increase which would require rid = 1/3 < 1. 

In the high wave numbers, depending on the value of 
the Schmidt number, two things can happen. Sc ^ 1 
means that buoyancy production cannot exist below k 
and cannot balance with viscosity to create gu (see 
FigH]). In this case, the velocity spectrum exponen- 
tially decreases after it goes through NLv. Otherwise, 
if 5*0 ^ 1 it is possible to prove that an equilibrium 
between buoyancy and viscosity exists by performing a 
linear analysis of pi2l3p with c{x, t) = cq + ci(fe, t) e~^^^ 
and u{x,t) — Uq + Ui{k,t) e""^^^ . This solution can be 
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thought of as the low wave number contributions of u 
and c with perturbation corresponding to higher wave 
number turbulent fluctuations. They can be plugged in 
equations (|H2l3p and after some algebra it is found that 



gj kjk 



fc2 



(11) 

c?(fc)e(*"°'=~'''^''*,(12) 



in the high Schmidt number limit. As a result and by 
definition of E{k) and Ec{k) it is then found 



E{k,t) 



ly^k^ 



Ec{k,t) 



(13) 



in agreement with Eq.Q. Therefore, the velocity spec- 
trum coincide with Eg^{k) in the high Schmidt number 
limit (see FigH]). 

As a conclusion, it can be stated that self similar- 
ity hypothesis together with equilibrium of spectral en- 
ergy transfer with buoyancy at low wave numbers con- 
strained velocity spectrum in a way incompatible with 
Kolmogorov mechanism. It does not provide the exact 
slope of the concentration spectrum Ud at low wave num- 
bers but assuming that it lies between 1 and 2 we obtain 
2 < ni < 2.5. Recent high resolution simulations, show 
that the velocity spectrum level remains close to the A 
line which would be in favor of rid ~ 1 ^md n/ w 2 
as a result. The most important result, here, is that 
velocity spectrum at low wave numbers is shown to be 
very sensitive to concentration spectrum. Nevertheless, 
a Kolmogorov mechanism is not ruled out at intermedi- 
ate wave numbers. It can be a transition process from 
low wave numbers to high wave numbers as in the BO 
mechanism. Most of the gravitational energy is injected 
through buoyancy at low wave numbers up to a critical 
wave number. Above it, injection of energy is negligi- 
ble and a Kolmogorov mechanism is possible. In this 
case, it does not involve that the turbulent dissipation 
e = 2i/ J dk k^E{k) should vary like e ~ As men- 
tioned in [m , this common belief that "the cascade rate 
is [...] equal to dissipation" in the self similar regime is 
always presupposed although it is a strong hypothesis for 
buoyant flows. It does not have to hold considering the 



importance of the buoyant production. It would not be 
in contracdiction with turbulent kinetic energy conserva- 
tion dt{u'^) = n — e. Indeed, since buoyancy production 
n ~ i and dt{u'^) ~ i, dissipation could vary slowlier 
without threatening energy conservation in the self sim- 
ilar regime. 
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